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One-shot decoupling is a powerful primitive in quantum information theory and was hypothesized 
to play a role in the black hole information paradox. We study black hole dynamics modeled by a 
trilinear Hamiltonian whose semiclassical limit gives rise to Hawking radiation. An explicit numerical 
calculation of the discretized path integral of the 5-matrix shows that decoupling is exact in the 
continuous limit, implying that quantum information is perfectly transferred from the black hole to 
radiation. A striking consequence of decoupling is the emergence of an output radiation entropy 
profile that follows Page’s prediction. We argue that information transfer and the emergence of Page 
curves is a robust feature of any multi-linear interaction Hamiltonian with a bounded spectrum. 


The explicit mechanism of black hole unitary evapora¬ 
tion is a major open problem in high-energy physics [1]. 
The evaporation of a macroscopic black hole is believed to 
be well described by a semiclassical result of Hawking [2], 
but it has also been argued that for a sufficiently old black 
hole the in-falling Hawking quanta (hereafter called the 
“partner modes”, as opposed to the radiated modes that 
constitute the Hawking radiation) must somehow resur¬ 
face in order to preserve the manifest unitarity of the 
black hole evolution. Page [3] argued for the latter using 
a description based entirely on calculating the entropy 
of the Hawking radiation as a subsystem of the joint 
pure state of a black hole and its radiation field, with¬ 
out proposing an explicit interaction that would enable 
this. In the modern language of quantum information 
theory, Page’s calculation implies that when the black 
hole entropy vanishes, it is decoupled from its surround¬ 
ings. Here we present a model in which black hole modes 
interact with radiation modes both inside and outside of 
the horizon in a fully unitary manner, and that allows us 
to explicitly calculate the entropy of the black hole over 
time as it gradually decouples from its surroundings. 

The modern approach to the resolution of the prob¬ 
lem started with Ref. [4]. These authors proposed an ab¬ 
stract mechanism for black hole evaporation based on the 
recently developed “decoupling” framework (see a gen¬ 
eral version in [5]). Decoupling theory is a fundamental 
technique to prove a plethora of quantum communica¬ 
tion protocols such as state-merging [6], the fully quan¬ 
tum Slepian-Wolf (FQSW) protocol [7] or the channel 
capacity coding theorem [8]. 

The central result of the one-shot decoupling theo¬ 
rem [9] implies the existence of a completely positive 
(CP) map Ta^b and a state pra where the following 


condition is satisfied 

\\Ta^b{UpraU^) — Db ® Pr\\i —> 0 (1) 

for almost all unitaries U. Here, ||A||i = Tr[(Af A) 1 / 2 ] is 
the trace norm, R is an auxiliary system, pr = Tr a[pra] 
and &b is a state specified later. For our purposes it will 
turn out to be sufficient to consider pra as a maximally 
entangled state (denoted ^/m). Then, decoupling Eq. (1) 
is equivalent to the conditions of the entanglement of the 
A subsystem of <&ra being transferred to E via another 
CP map Ta^tE called a complementary channel to Ta->b ■ 

In what follows we examine decoupling induced by a 
boson trilinear operator describing the interaction of a 
quantized black hole’s modes with those of radiation in¬ 
side and outside of the event horizon, in the special case 
when Ta^e is an identity, something that, crucially, is 
not an assumption of the model but rather an unexpected 
and natural consequence of the interaction Hamiltonian 
we use. 

In the semiclassical description, the outgoing Hawking 
radiation shrinks the black hole but the exact mecha¬ 
nism of how the black hole loses its mass over time is 
beyond the scope of Hawking’s original derivation, as in 
that approach states from past infinity are mapped to 
future infinity via a Bogoliubov transformation. In order 
to restore an explicit time-dependence to the evaporation 
process, we can use the by now well-known insight [10, 11] 
(see also [12, 13]) that any solution to a time-dependent 
quantum harmonic oscillator can be understood in terms 
of a Bogoliubov transformation that connects the initial 
to the final frequency, and in particular can be imple¬ 
mented via a time-dependent Hamiltonian of the “squeez¬ 
ing” type in which the coupling strength r{t) plays the 
role of the classical time-dependent driving field (we set 
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n = i) 

H = — be) . (2) 

Here, the annihilation operators b and c annihilate ra¬ 
diation modes within and outside of the event (or ap¬ 
parent) horizon, respectively. Such a description can 
be used to derive Hawking radiation from a collapsing 
shell of matter, for example [14-16]. However, in such a 
time-dependent semi-classical description, the black hole 
degrees of freedom are not explicitly quantized and the 
backreaction of the radiation on the black hole cannot be 
followed. 

The dynamic model for the evaporation of a black hole 
that we analyze here can be seen as an an extension of 
Hawking’s semi-classical approach [17] where the black 
hole internal degrees of freedom are explicitly quantized 
to allow us to track the black hole quantum entropy 
over time. This model is explicitly ad hoc: the quan¬ 
tum interaction is not derived from a quantum field the¬ 
ory that couples gravitational modes to radiation modes. 
Rather, it is an extension of the semi-classical Hamilto¬ 
nian (2), and reduces to it when the quantum black hole 
modes are replaced by their classical expectation values. 
In the static infinite-time limit, the quantum Hamilto¬ 
nian (shown below) therefore implements the standard 
Bogoliubov transformation between past and future in¬ 
finity, and reproduces the celebrated Hawking results as 
has been shown previously [11, 18, 19]. 

The Hamiltonian we study is the so-called boson tri- 
linear Hamiltonian 

Htri = *r(afeV — a^bc) (3) 

often encountered in quantum optics and condensed mat¬ 
ter physics [20-22] or quantum optomechanics [23]. In 
quantum optics, the operators b and c refer to the signal 
and idler modes, respectively, while a refers to the pump 
mode. The intriguing formal analogy between black hole 
radiance and parametric amplification has been noted be¬ 
fore [18, 19, 24]. As a model of black hole and radiation 
interactions, the operator a in (3) annihilates black hole 
modes, whereas c and b are defined as in (2). 

The black hole mass is a function of the input state 
n a = (a'a) as well as of the coupling strength r (called 
the gain or squeezing parameter in quantum optics). In 
the semi-classical limit, the relationship between r(t) in 
(2) and black hole parameters is simply (in units where 
h = c = G = k= 1) [18] 


21n[tanh 1 (r(t)t)] 

where T(t) is the time-dependent Hawking temperature 
T(t) = (87rM(t)) _1 and M(t) is the black hole mass. 

The initial black hole state can be pure or mixed. For 
simplicity, we restrict our analysis here to the pure state 
scenario and a black hole in initial state | n) a (n modes), 


but our results naturally carry over to an arbitrary super¬ 
position and any mixed initial black hole state with cor¬ 
responding nonzero initial entropy as we discuss below. 
We now define the time-evolution operator on a quantum 
state l'k(ti)) as |4'(t 2 )) = U(t 2 ,ti) |4/(ti)) where 

t/(t 2 ,fi) = Te- i ft H{t ' )dt '. (5) 

The symbol T stands for Dyson’s time-ordering oper¬ 
ator and H(t ) is the Hamiltonian describing the uni¬ 
tary evolution of the quantum state. We can then 
write the time-evolution of the black hole initial state 
|*(0)> = |n) o |0) 6 |0) c = |n) a |0) 6c as 

!*(*)> = U(t,0) l’F(O)} = Te-^o ln)a | 0)fec) (6) 

where H tr] is the trilinear Hamiltonian (written as a sum 
over modes) 

OO 

Htri= E ir ^k( t )( a kblcl-alb k c k ). (7) 

k =—oo 

In the following, we will focus on a single mode k with 
energy cu*,, and omit the index k for convenience. 

To calculate (6), we introduce small time slices At and 
discretize the path integral so that with t = NAt 

N 

U(t ) = T e~’ Jo n^'Ut' » H e~ iAtHi , (8) 

i= 1 

where the *-th time-slice Hamiltonian Hi = «ro(a6|cj — 
a^biCi) acts on the black hole state and the i-tli slice of 
the be Hilbert space |0) b . c .. The initial value ro sets the 
energy scale, and we simply set ro = 1 in the following. 
We now apply (8), so that 

N 

mt)U c = W | n) a |0> bc = n e~ iAtH > I n) a |0> bc , (9) 

i =1 

where |0) bc = |0) bjvcjv ® ... <8> |0) 6ici , = e ~ iAtHi and 

W = IfW o ... o The first time slice is 

n 

l*(l)> = W^\n) a |0) blCl = £ \ j) a \n j) bici 

3=0 

(10) 

while l’F(t)) is given by 
\m)a’ bc = E U nJ---U%l ljN 

31---3N 

X \jN)a\jN-l - jN)b N c N ■ ■ ■ \n ~ jl)b lCl - (11) 

Note that approximating the path integral using a sin¬ 
gle time-slice—the static path approximation (SPA), see 
e.g. [25, 26]—can yield good results at very low tem¬ 
peratures, but disregards self-consistent temporal fluctu¬ 
ations. Indeed, SPA calculations of the black hole en¬ 
tropy using the trilinear Hamiltonian lead to an oscillat¬ 
ing behavior of the black hole entropy [18, 19]. The black 
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FIG. 1. The action of the discretized path integral W *> is 
sketched. The braces indicate the entropy calculation of the 
radiated modes at each time-slice. 


hole von Neumann entropy can be calculated using the 
marginal density matrix •d a '(t) = Tr bc [|^ , (t))(T'(t)| a 'bc], 
which can be written entirely in terms of the coefficients 
\Uij\ 2 , which describe the probability to observe j modes 
of the black hole given that in the previous iteration we 

n 

found i modes) $ a / (t) = Pjn I Jjv )<JJV | where 

iw=o 


sity matrix $„/ (for N 0) is extremely close to a pure 
vacuum state |0) a , for all examined input basis states 
| n). Note that the case n = 50 gives rise to an extremely 
large Hilbert space, and to obtain a reliable Taylor ex¬ 
pansion for At = 1/15 would require expansion to order 
up to 500. A straightforward numerical evaluation of 
the 2 500 terms is of course impossible, but the method 
developed in [? ] makes the calculation tractable on a 
High-Performance Computing Cluster. We have also per¬ 
formed a calculation of the integral where partner modes 
are allowed to interact with the black hole a second time, 
breaking the symmetry between partner and Hawking 
modes. Such a calculation is far more complex, but does 
not change the shape of the Page curves appreciably. 

It is important to stress that even though we studied 
the effect of W only on basis states |n), we can extend 
this conclusion to an arbitrary input pure state ipa = 
J2k=oPk |fc) a - It is sufficient to show decoupling of the 
spanning basis set. Thus, ultimately the action of W can 
be reformulated as inducing 

|fc)„|0) 6c ^|^)U c «|0) a ,®^ c (|fc», (13) 


V = V \U {1) \ 2 I U {N - 1] | 2 | U (N) I 2 

PjN / ^ \ U nj 1 \ * • * \ U 3N-23N-1 1 \ U jN-13N I * 

( 12 ) 

To follow the time evolution of the black hole entropy 
H n (t) = — Tr [$o'(t)log$a'(t)] given the initial state of 
the black hole is a pure state with n quanta, we need 
to evaluate the n coefficients pj N above. This is a dif¬ 
ficult calculation because unlike the Hawking isometry, 
the unitary operator is not known to be factorizable 
in a simple way. The usual formal factorization formu¬ 
las [27-29] are infinite operator products and often are 
not suitable for practical calculations. A common route 
to obtain the time evolution is to expand the exponential 
around r$t ss 0 but this allows only a short-time analy¬ 
sis of the evolution operator. For greater r 0 t, the Taylor 
expansion is prohibitively inefficient. However, a method 
developed in Ref. [? ] enables us to find the action of the 
unitary operators even for large values of rot. The action 
of WW acquires an interesting combinatorial interpreta¬ 
tion in terms of an integer lattice known as a generalized 
Dyck path [30] if it acts on any state generated by the 
repeated action of atfc' on a ground state |gr), defined 
by a)bc |gr) = 0. The basis elements {\k 00) abc }£ =0 span¬ 
ning our input Hilbert space of are all ground states 
of Hi. 

We evaluated the output entropy for black holes with 
initial quanta n for 0 < n < 50 using the discretized path 
integral (11). The resulting entropy as a function of the 
number of time slices (at constant small At) is shown in 
Fig. 2 for n = 5, 20 and 50. An entropy curve strikingly 
similar to the one predicted by Page [3] emerges in the 
limit of large t (several thousand time slices). This shape 
is consequence of the fact that the final black hole den¬ 


where 0 < k < n and $b c (|fc)) = Tr a / [|\E>(f))('F(f)| a / 6c ]. 
Hence $b c ( \k)) is necessarily pure and mutually orthogo¬ 
nal for all k. Since |0) a is a “constant” state (independent 
of k), it immediately follows that an arbitrary superpo¬ 
sition ip a will be transferred to the Hawking modes and 
their partners (bc). These are precisely the circumstances 
captured by the decoupling theorem in Eq. (1) upon iden¬ 
tifying A ±=, a, E i=f bc and B ±=; a'. The presence of an 
input bipartite state Qra instead of if )a is just a conve¬ 
nient reformulation so that if the decoupling (transfer) 
of is successful, it can be used to perfectly teleport 
any ^ A - 

As is implied by (13), the black hole evolution ulti¬ 
mately turns out to be an erasure map T a ^a' ■ 4>a —>• 
Tr[V’a]|0)(0| a /, which we interpret as complete evapora¬ 
tion. This immediately implies that the map T a ^>bc must 
be the identity as announced earlier. The decoupling 
theorem then implies [5] that the input state (and so the 
quantum information it carries) has been nearly perfectly 
transferred to the bc Hilbert space. What we found there¬ 
fore is an explicit (and simple) realization of the one-shot 
state merging protocol [5, 31] (generalizing Ref. [6] stud¬ 
ied in the asymptotic limit) and thus the action of T a ^,a' 
in (1) turns out to be a simple case of the FQSW pro¬ 
tocol [7] (the erasure map is essentially the trace map 
which is a special case of the partial trace appearing in 
FQSW). 

The question of how exactly the information is encoded 
in the outgoing radiation described by the density matrix 
i9bc is non-trivial. The mapping T a ^bc is a highly convo¬ 
luted identity, something that is expected from the de¬ 
coupling mechanism. Hence the information is not read¬ 
ily available for a read-out—the state occupying the bc 
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FIG. 2. Black hole entropy H n (t) = H[$ a i] = H[&b c (t )] as a 
function of the iteration number N for initial states | n) |0) b 
with n = 5 (orange), n = 20 (red) and n = 50 (pink), with 
t = AtN. Logarithms are base n + 1 so that each black hole 
entropy H n (t ) < 1. We use ro = 1 and At = 1/15. The curves 
are qualitatively unchanged if At = 1/25 and using a com- 
mensurately larger number of time slices (data not shown). 

Hilbert space is scrambled, or, in quantum information 
theory jargon, it is a highly entangled quantum code. 

The output state is highly entangled and so an ob¬ 
server having access only to the b or c mode cannot learn 
anything. This is a direct consequence of the symmetry 
between the modes b and c in the boson trilinear Hamil¬ 
tonian given by Eq. (3), also seen in the action of the 
unitary W W and W in Eqs. (10) and (11), respectively. 
This is a textbook example of how the impossibility of 
perfect cloning prevents an observer from obtaining any 
information from i or d c alone. If there existed a com¬ 
pletely positive map that recovers i\) a from db only, the 
same operation could recover it from d c . But that would 
result in two nearly perfect copies of ijj a which is pro¬ 
hibited by quantum mechanics. Using the language of 
quantum Shannon theory, the quantum capacity of the 
channel A4 a _>b(c) must vanish [32]. 

The trilinear Hamiltonian Eq. (3) that allowed us to 
carry through this analysis is, to put it bluntly, a guess. It 
does not follow from an interacting quantum field theory 
of gravity (such as a locally-flat conformal field theory 
coupled to a dilaton [33] (see [34] for additional refer¬ 
ences). However, it is manifestly an elegant extension of 
the semi-classical approach, and it is difficult to imagine 
that a consistent renormalizable quantum field theory of 
gravity would not, in an effective limit, reduce to pre¬ 
cisely an interaction as we describe. In fact, the mech¬ 
anism of information conservation implied by our model 
is an example of a much larger class of more compli¬ 
cated Hamiltonians with the proper semiclassical limit 
V Ui |vac) = 1/coshr Wi J2m=o tanh m |m) b . |m) c .. In¬ 
deed, it may be possible to use decoupling as a con¬ 
straint to impose conditions on the form of the Bogoli- 
ubov transformation that connects the in-vacuum with 


the out-vacuum. It is clear from general arguments [35] 
that at a minimum extra degrees of freedom are needed 
that are coupled to the Hawking modes, and in our model 
these degrees of freedom are the quantized black hole 
modes a. 

To illustrate the mechanism concretely, we set = 
|3) a and observe the behavior of the probabilities given 
by Eq. (12) in Fig. 3. The black hole state slowly “drifts” 



FIG. 3. The schematic evolution of the probability Eq. (12) 
of the black hole density matrix t? a / (t) for the input state |3) a . 
Time goes from left, the boxes are Fock bases (|0),..., |3), top 
to bottom) and the lines represent the contributions from the 
previous iterations. The vacuum gets the majority of contri¬ 
butions in every iteration (three steps shown here) capturing 
the black hole evaporation. 

towards the vacuum state (the top line) even if the tran¬ 
sition amplitudes for the highest Fock state are close to 
one, that is, the least favorable conditions for the tran¬ 
sition to the vacuum state. This is a direct consequence 
of the spectrum boundedness of the trilinear Hamilto¬ 
nian Eq. (3) - the number of bosons n in the a mode 
that W = exp [— iH tl i\ can “distribute” to the be modes 
is naturally equal to or lower than n. Crucially, any 
spectrum-bounded and possibly multilinear Hamiltonian 
displays exactly the same behavior, and only the infor¬ 
mation rate differs. As a consequence, the information 
transfer mechanism-and therefore the Page curve-are ex¬ 
pected to remain virtually unchanged even for much more 
complicated internal dynamics. In other words, the Page 
curve is robust and universal. 

The general features of information transfer between 
quantum subsystems as described by Page [3] (see also 
more recently [4, 36, 37]) have been discussed intently in 
the past. We believe our calculation is the first explicit 
realization of these previous abstract arguments, point¬ 
ing towards decoupling as a mechanism for the resolution 
of the black hole information loss problem. Finally, we 
would like to point out that the explicitly time-dependent 
black hole evaporation process we described throws new 
light on the discussion of putative firewalls [38, 39] for old 
black holes. The usual line of reasoning for firewalls based 
on the monogamy property of entanglement does not ap¬ 
ply here since the black hole state | v E'(t)) a , bc from (11) is 
tripartite entangled. 

CA thanks Steve Giddings for discussions. This work 
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